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FINITE AMPLITUDE WAVES IN FLUID-FILLED ELASTIC TUBES: 


WAVE DISTORTION, SHOCK WAVES, AND KOROTKOFF SOUNDS 
By Richard M. Beam 
Ames Research Center 


SUMMARY 


The Lagrangian form of the one -dimensional equations of motion for a 
fluid-filled elastic tube is developed. The resulting nonlinear equations are 
combined -with a general nonlinear internal pressure cross-sectional area rela- 
tion for the tube to obtain the finite amplitude "wave equation." After 
determination of the tube pressure-area relation which leads to a linear (non- 
distorting) wave equation, a solution for the general (distorting) wave equa- 
tion is presented. The solution for the general wave equation is applicable 
for simple (nonreflecting) waves in tubes with general pressure-area relations 
The solution is then used to develop criteria for the steepening and 
nonsteepening of finite amplitude waves . 

The theory is further developed to include the calculation of wave distor 
tion. The criterion for a sharp wave is derived and the critical length of 
tube required for a sharp wave to develop is determined. The propagation of 
the resulting sharp wave or shock wave is considered. The velocity and the 
energy loss formulas for the shock wave are developed. 

Finally, the analysis and experimental results are used to support new 
hypotheses for the mechanism and source of the Korotkoff sounds which have a 
major role in the indirect determination of blood pressure. A limited amount 
of experimental data (from Brarrrwell) is presented to correlate quantitatively 
the analysis and the new hypotheses. 


INTRODUCTION 


The study of the dynamics of fluid-filled elastic tubes, or cylinders, 
has received extensive interest over the past two centuries. The greatest 
bulk of papers and books on the subject has appeared within the last 10 years 
in the fields of biomechanics 1 and aeromechanics (see the bibliography of 
ref. 2). Studies in aeromechanics are, in general, concerned with the inter- 
action of liquid propellants and their elastic containers while the studies in 
biomechanics are concerned with the interaction of blood and arteries or veins 
Although the physical analogies in this report refer to biomechanics problems, 
the application of the theory is not restricted to one field. 


1 An excellent historical and theoretical review including an extensive 
bibliography was given by Shalak in 19 66 , reference 1. 


For mathematical simplicity, most investigators have restricted their 
analyses to linear theory. The solutions to the nonlinear equations have been 
limited either to numerical integration utilizing the method of characteris- 
tics (ref. 3) or to the perturbation methods (ref. l) . Although both methods 
of solution provide useful information, neither provides the general 
information available from a closed-form solution. 

The purpose of this report is twofold; first , to present an analysis for 
the propagation of a simple (no reflection) wave of finite amplitude in a 
fluid-filled elastic tube. Of particular interest are the distortion of a 
smooth wave until a sharp wave develops and the propagation of the resulting 
shock wave. Although several writers have speculated on the development of 
shock waves (refs. 1, 4, and 5) in fluid-filled elastic tubes, this is, to the 
knowledge of this author, the first presentation of a theory for the 
development and propagation of the shock waves. The experimental verification 
remains to be done. 

The second purpose of this report is to present new hypotheses for the 
mechanism and source of the Korotkoff sounds which were the motivation for 
this study. The hypotheses concern the source of the Korotkoff sounds that 
occur when blood pressure is measured with a pressure cuff. Briefly, the 
mechanism hypothesis states that the sounds result from the development of a 
sharp, or possibly, a shock pulse wave in the compressed portion of the artery 
under the inflated cuff. The source hypothesis states that the source of the 
"sounds" is the response of the "listening" system to a sharp pressure wave. 


SYMBOLS 


A(x,t) 

A e 

Ai,A2 

C 

c(A) 


Cl 





internal cross-sectional area of tube 

equilibrium internal cross-sectional area of tube 

discrete values of cross-sectional area of tube on each side of 
shock wave 

wave velocity of finite amplitude nondistorting wave 

spatial wave velocity as function of area A (the argument may be 
other than A) 

spatial velocity of shock wave 

shock -wave velocity in compressed coordinate system 

a constant (eq. (59)) 

energy loss coefficient for shock wave 

function defined by equation (27) 
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f(a) 

fz(A) 



h [ t ‘^) 
L c (p*) 

^cm 

p(x,t) 

Pe 

P* 

P2>Pl 

t 

"to 

u(x,t) 

v(x,t) 

V 2 ^Vl 

X 

J 

<x(x,t) 

e(x^t) 


€l> e 2. 


pressure cross-sectional area relation for tube 

pressure cross-sectional area relation for tube which leads 
to a linear wave equation 

function defined by equation (25) 


pressure spatial time relation for simple wave 

length of tube required for a sharp wave to develop at 
pressure p* 

minimum length of tube required for a sharp wave to develop 

pressure differential across tube wall 

(Subscripts i and o denote inside and outside tube 
P = Pi - Po- ) 

differential pressure across tube wall at equilibrium 
condition 

(Subscripts ie and oe denote inside and outside tube 
pressure p e = p ie - P oe -) 

pressure at which pressure gradient becomes infinite 

discrete pressures on each side of shock (px is also used as 
constant in eqs. (58) through (63)) 

time parameter 

discrete time at which pressure equals p* (eq. (52)) 

displacement of particle which had coordinate x at time 
t = 0 

particle velocity (eq. (3*0) 

discrete particle velocities on each side of shock 

spatial coordinate (fig. l) 

spatial coordinate (fig. 3) 

dimensionless relation for tube cross-sectional area 
(eq. (32)) 

dimensionless relation for tube cross-sectional area 
(eq. (20)) 

discrete values of € on each side of shock 
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\i particle -to -particle propagation rate (eg.. (36)) 

p fluid density 

w frequency constant 

(“) dummy variable of integration 


I. LAGRANGIAN EQUATIONS OF MOTION 


Derivation of Equations of Motion for Large Displacements 

The purpose of this section is (l) to develop the large displacement equa- 
tions of motion (Lagrangian form) for a fluid-filled elastic tube, ( 2 ) to 
establish the conditions for a nondistorting propagating wave, and (3) to 
obtain the finite amplitude wave propagation properties. In order to simplify 
the mathematics and to emphasize the physical phenomena relevant to this 
report, some gross assumptions regarding the physical system have been made, 
namely: 

(a) The fluid is incompressible and inviscid. 

(b) The fluid flow is one -dimensional and the radial inertia is 
negligible . 

(c) The mass of the tube wall is negligible (i.e., the kinetic energy of 
the wall is small when compared to the kinetic energy of the fluid). 

(d) The flexural rigidity and prestress of the tube in the axial direc- 
tion do not affect the stiffness of the tube. This is equivalent to the 
assumption that each element along the length of the tube behaves as an 
isolated element disconnected from its neighbors. 

(e) Under equilibrium conditions the tube has uniform (not necessarily 
circular) cross section along the length of the tube. The elastic properties 
are also uniform along the tube. 

The validity of these assumptions in the analysis of Korotkoff sounds 
will be reviewed in section IV of this report. 

Euler (ref. 6 ) developed two forms for the equations of motion of a fluid. 
One considers the velocity, pressure, and density of all points in the space 
occupied by the fluid for all time, while the other considers the history of 
every element of the fluid for all time. The two forms are generally desig- 
nated as the Eulerian and Lagrangian forms, respectively. For most fluid flow 
problems, the Eulerian form is more tractable than the Lagrangian, and the 
Lagrangian form becomes quite awkward for two- and three-dimensional flow. 

For some special one -dimensional flows, however, the Lagrangian form has some 
advantages which will be evident in the following analysis. 
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Consider the equilibrium of an element of fluid which had original length 
5x. Let u(x,t) be the displacement of the particle which had coordinate x 
at time t = 0 (see fig- l) - The x coordinate system is then an inertial 




8x 



coordinate system along the axis of the tube. In addition, define 

A(x,t) cross-sectional area at time t of the fluid element which had 

coordinate x at t=0 (i.e., the cross-sectional area of the 
element at x + u) 

A e original (t = 0) cross-sectional area of the fluid element 

p(x,t) pressure at time t in the fluid element which had coordinate x 

at t = 0 (i.e., the pressure in the element at x + u) 

p fluid density (a constant for incompressible fluid) 

The force on the left side of the deformed element is 


p(x,t)A(x,t) 


( 1 ) 


while the force on the right side is 

-p( x -l- 5x,t)A(x + Sx,t) (2) 

In addition, the forces acting on the sides of the element have an axial 
component 


p( x, t ) + p( x + 5x,t) 


[A(x,t) - A(x + 5x,t)] 


(3) 


The forces acting on the surfaces of the element (eqs. (l), (2), and (3)) must 
be offset by the inertial force of the accelerating element 


pSxA e 


![ 


5 2 u(x,t) <3 2 u(x + Sx,t) 






(4) 
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Note that expressions (3) and (4) are only approximate for finite Sx hut are 
exact in the limit as Sx approaches zero. Equating the external forces to 
the inertial forces acting on the element one obtains 


p(x,t)A(x,t) - p(x+ 5x,t)A(x + Sx,t) — - - ^4 ~ - [A(x,t) - A(x+Sx,t)] 


= pSxA e - 


S 2 u(x,t) 5 2 u(x + 5x,t) 

St 2 St 2 


(5) 


If the terms of equation (5) with argument x + Sx are expanded about x 
p(x,t)A(x,t) 


p(x,t) + Sx — - + 0(Sx 2 ) 


A(x,t) + Sx — + 0(Sx 2 ) 


Sx 


1 

2 


p(x,t) + p(x,t) + Sx ^4 ~ — “ + 0(8x 2 ) 

ox 


A(x,t) - A(x,t) - Sx - 0(Sx 2 ) 


p5xA e ^ 


'S 2 u(x,t) S 2 u(x,t) 

St 2 + St 2 + 


Sx 


5 2 5u(x,t) 
at 2 Sx 


+ 0(Sx 2 ) 


( 6 ) 


where 0(Sx 2 ) represents terms involving terms of second and higher order in 
Sx and, after terms of 0(Sx 2 ) are regrouped, equation (6) may be written 

- | d P d^ t ^ A(x,t)Sx + 0(5x 2 ) = p5xA e ^ ^* 2 ^ (?) 

Since Sx may be taken as small as one chooses, in the limit as Sx -*• 0 
equation (7) becomes 


dp(xjt) 

Sx 


A(x,t) 


pA e 


S g u(x,t) 

St 2 


(8) 


In addition to equilibrium of the deformed element, the mass, or in this 
analysis (incompressible fluid) the volume, must be conserved. Equating the 
volume of the undeformed element to the volume of the deformed element 
produces 

5xA e = [Sx + u(x + 5x,t) - u(x,t)] ^ [A(x,t) + A(x + Sx,t)] (9) 

which is approximate for finite Sx but exact in the limit as Sx approaches 
zero. 
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Pie pressure inside the tube when tube is in 
equilibrium with outside pressure p oe 
and has inside cross-sectional area Ae • 


After expansion of equation (9); 
the limit as Sx -* 0 is 


A e - 



5u(x,t ) 
dx 


A(x,t) 


( 10 ) 


and substitution from the mass conserva- 
tion equation (10) into the equilibrium 
equation (8) leads to an equation relat- 
ing the pressure and the displacement 


dp(x,t) 

“ ^ " P 


du(x,t) 

Sx 


5 2 u(x,t ) 
St 2 


( 11 ) 


The preceding equations are inde- 
pendent of the properties of the tube 
wall. The following analysis will be 
restricted to the special case where 
the internal cross-sectional area of 
the tube and the pressure difference 
across the tube wall are related by a 
single- valued function f, in particular 


p oe pressure outside the tube when tube is in 
equilibrium. 

p e differential pressure (p^ - p Q ) across the 

tube wall when tube is in equilibrium 
with internal cross-sectional area Ae • 

A e equilibrium internal cross-sectional area 

of tube . 

p change in differential pressure across 

tube wall from p g (subscripts i and o 
refer to inside and outside tube) . 

A cross-sectional area of the tube wall 

corresponding to pressure p. 


Figure 2.- Pressure-area relation for typical 
elastic tube. Dashed axes indicate pressure- 
area relation applicable for perturbation 
about equilibrium conditions (p e >A e ). 


P = f(A) (12) 

The qualitative shape for a typi- 
cal elastic tube pressure -area curve is 
illustrated in figure 2. Note that p 
of figure 2 represents the differential 
pressure across the tube wall but the 
p(x,t) in equations (l) through (ll) is 
the pressure in the fluid. This appar- 
ent inconsistency in notation is recti- 
fied by letting the outside pressure 
change be zero (p 0 = constant) through- 
out the remaining analysis. Note also 
that the pressure-area curve depends on 
the initial cross-sectional area, A e , 
and pressure, p e . 


The pres sure -area curve of fig- 
ure 2 includes the full range of cross-sectional area. With large negative 
values of p 0 (pj_ e « p Q e) the tube is completely collapsed (A e = 0) . As the 
external pressure is decreased (or internal pressure is increased) the area 
gradually increases until the pressure difference approaches the classical 
buckling pressure for the tube (ref. 7)- In practice,, the tubes do not 
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exhibit instability characteristics since they are not perfect cylinders. 2 
Rather, they exhibit regions of large change in area for small changes in 
pressure. As p e becomes positive the cross section of the tube becomes 
quite circular and remains so until the pressure becomes sufficiently high at 
which time the tube reaches an unstable equilibrium configuration. The pres- 
sures at which this later type of instability occurs are outside the range of 
interest of this report. 

With the pressure-area relation defined by equation (12) the derivation 
of the equation for u(x,t) may be completed. Since the pressure is known as 
a function of A then 

dp(x,t) df (A) dA(x,t) , r / N-, dA(x,t) /, v 

Sx “ cLA dx LMX,t;J 8x C13 ‘ 

but from differentiation of both sides of equation (10) with respect to x 


A Jj,( x ; t) 

dA(x,t) e 8 X 2 



or from (13) and (l4) 


dp(x,t) 

8x 


-f *[A(x,t)]A 6 


1 + 


8 2 u(x,t) 

5x g 

6u(x,t) 

8x 


Equation of 


6p(x,t) 

8x 


from equations (ll) and (15) leads to 


5 2 u(x,t) 


f *[A(x,t)]A e 


1 + 


6x 2 

5u(x,t) 

8x 


5 u(x,t) 
St 2 


(14) 


( 15 ) 


(16) 


Equation (l6) is the general nonlinear wave equation for large- amplitude waves 
and a general nonlinear tube pressure-area relation. In order to obtain an 
equation for u(x,t) alone, A can be eliminated from equation (l6) by 


2 Although the collapse of an artery differs from that of a rubber tube 
(see ref. 8), the essential features of figure 2 are the same for both. 
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substitution of A from equation (10) . The resulting equation will, in 
general, be nonlinear (the complexity depending on the form of f*(A)) but 
susceptible to analysis. 


Condition for Nondistortion of Large -Amplitude Wave 


It will prove advantageous to eliminate the bracketed term in equa- 
tion (l6) by substitution from equation (10) 


71“ f , [A(x,t)]A 3 (x,t) 
A e P 




(17) 


Note that equation (17) has the form of the linear wave equation if 


f ^[A(x,t) ]A 3 (x,t) = constant 


and that the corresponding wave velocity is C provided 


~h~ fj[A(x,t)]A 3 (x,t) = C 2 ( 18 ) 

A e P L 

The required f(A), denoted fj(A), is found from equation (l8), in the form, 


df z (A) = C 2 pAg |§ 


therefore 


fl(A) = C 2 pA| + const 

The constant of integration is obtained from fj(A e ) = 0, therefore, 

f Z (A) = ^ 


, A e V 

1 - ( f 


(19) 


If the pressure -area relation has the particular form prescribed by equa- 
tion ( 19 ) > the linear wave equation will prevail and simple waves with finite 
amplitude will propagate without distortion and with velocity C. It is inter- 
esting to note that a longitudinally constrained rubber tube with circular 
cross section has precisely the pressure-area relation of equation ( 19 ) 

(finite amplitude displacements with cross section remaining circular). 

The nondistortion of large -amplitude waves in rubber tubes with circular 
cross section was pointed out previously by Olsen and Shapiro (ref. 3)« Using 
the Euler ian form of equations, they introduced the pres sure- area curve far a 
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rubber tube (circular cross section). After developing the equations for the 
characteristic curves, they considered a finite amplitude wave traveling in 
one direction and with the appropriate argument concluded that "every part of 
the wave moves with the same speed and that the wave neither steepens or 
flattens. This extraordinary result depends on the particular pressure -area 
relation for the tube." 

The above analysis (eqs. (17) through (19)) proves that the only form of 
pres sure -area curve which yields a linear wave equation (nondistorting waves) 
for finite amplitude waves is given by equation ( 19 ). 


Finite Amplitude Wave Distortion 

The primary concern of the present investigation is the distortion of 
finite amplitude waves as they propagate. Return now to the general nonlinear 
equation of motion (eq. (l6)). From equation (10) 


where 


then 


but 


therefore 


e(x,t) 


A e 

A(x,t) "*■ 


or 


A(x,t) 


A e 

1 + e(x,t) 


e(x,t) 


du(x,t) 

dx 


f'(A) 


df(A) 

dA 


df ( £ ) M = f »/-\ de 
de dA K J dA 


de ^e (1 + e) 2 
dA ~ A 2 A e 


f’(A) 


-f’(e) 


(1 ± e ) 2 

A e 


( 20 ) 


( 21 ) 


( 22 ) 


Substitution of f*(A) from equation (22) into the general equation of 
motion (l 6) leads to 


(1 + e) 


ci 2 u( X,t) 5 5 u(x,t ) 


at' 


at' 


(23) 


where e is, of course, a function of x and t defined by equation (2l). 
Equation (23) could, of course, be written 
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(23a) 


or, more generally. 


f Su(x, 

L s x 


t) 


p 1 + 


Su(x,t) 

Sx 


I S 2 u(x,t) _ S 2 u(x,t) 
St 2 " ^fc 2 


where 


g 


g 


Su(x,t) 

S 2 u(x,t) 

. S 2 u(x, 

t) 

Sx 

Sx 2 

£fc 2 


-p t 

Su(x,t)l 


Su(x,t) - 

I 

. Sx ^ 


- Sx 

“ " . r. 

. Su(x.t)l 


1 + 


dx 


(24) 


(25) 


Consider the propagation of a wave to the right (x positive) into the 
undisturbed tube. The solution must satisfy the differential equation (24) in 
addition to the conditions of the undisturbed tube ahead of the wave, namely 


du(x,t) 

dx 


= e(x,t) = 0 , 


Su(x,t) 

St 


( 26 ) 


In accordance with the development of Earnshaw (refs. 6 and 9 ) i we seek 
solutions of the form 


Su 

St 



(27) 


where F is a function to he determined. One may deduce from differentiation 
of equation ( 27 ) 


S /Su\ _ _s_ 
bt \§t J St 


F 



F* 


/ou\ S /ou\ 
\bxj St \Sx/ 


F* 


Su\ _S_ /SuN 
S xj Sx \Sty 



Su\ 2 S 2 !! 
Sx/_ Sx 2 


Therefore 


S 2 u 

St 2 



/'SuVl 2 S g u 

WJ Sx 2 


Equation (27) provides a solution to equation (24) if 


( 28 ) 
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( 29 ) 



Therefore equation ( 29 ) is the equation for f (M) ° rF(e) - 


Integration of equation ( 29 ) leads to 



However, if the integration is begun in the undisturbed region, one obtains 
from the conditions ( 2 6 ) and equation ( 27 ) 




61 = ^ 


therefore 


F(e) = ± 



de 


( 30 ) 


The particle velocity as a function of area can thus he obtained with the aid 
of equations (20), (22), ( 25 ), and (30) 



A compression wave (A > A e ) traveling to the right (x positive) produces 
positive particle velocity and a decompression wave (A < A e ) traveling to the 
right produces negative particle velocity. The converse is true for waves 
traveling to the left (A > A e > du/dt < 0; A < A e > du/dt > 0) ; therefore the 
positive sign preceding the integral of equation (3l) corresponds to waves 
traveling to the right and the negative sign corresponds to waves traveling to 
the left. 


For future reference in the evaluation of the wave distortion, it will be 
useful to introduce yet another parameter a defined as 
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a 


(32) 


1 

2 


1 



- | (26 + e 2 ) 


Note that a has "been chosen such that the linear wave equation prevails if 
the pressure is a linear function of a (i.e., if f(a) ~ a, see eq. ( 19 )). 
The particle velocity as a function of a is obtained from equations (30) , 
( 31 ), and ( 32 ) 


Su(x,t) 

St 



f(e) ~ 
p(l + e) 


1/2 


de = + 



+ 


p(l - 2a) 


1/ 2 


da 


( 33 ) 


The sign of the integral is, in accordance with the discussion of 
equation (3l), positive for waves traveling to the right and negative for 
waves traveling to the left. 

The final objective of this section is to evaluate the spatial velocity 
of propagation of a disturbance and the distortion of the disturbance. The 
procedure s will be to evaluate the rate of propagation from particle to parti- 
cle and with this result to consider the absolute velocity of propagation, or 
the 'wave velocity, of the disturbance. 

Let v(x,t) denote the particle velocity, that is, 


v(x,t) 


Su(x,t ) 


(34) 


The velocity of a particle at x at time t will be the same as the velocity 
of the particle at x + Sx at time t + St^ provided 


v(x,t) = v(x + Sx, t + St) 


or 


Msili Sx + Miiii , 0 

ox St 


but from equations ( 27 ) and (34) 

Su(x,t) 

therefore 


St 


= v(x,t) = F 


Substitution of 


S 2 u 


St 


/Su\ S u(x,t) „ s u( X, t ) 
\fr) Sx + ^ 

2 from equation (28) produces 


= 0 


(35) 
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F f 



d 2 u(x,t ) 




Sx + F 1 



which is true for all values of x and t if 



= 0 


Therefore if the rate of propagation from particle to particle 
velocity v is denoted |i , then 



of the 




Sx 

St 



(36) 



Figure 3.- Definition of the spatial coordinate 


Let the absolute location of a 
particle with velocity v(x,t) be 
y(x,t) (see fig. 3)> that is, 

y(x,t) = x + u(x,t) (37) 

then the absolute location of a parti - 
cle with velocity v at time t + St 
y.will be 


y(x + pSt, t + St) 


y(x,t) + ^ P&t + 


6y 

dx 


St 


Thus the absolute spatial velocity of particle velocity v (denoted 


Dy/Dt) is 


Dy _ 6y dy 
Dt dx 1-1 dt 


(38) 


Note that this represents the propagation rate of particle velocity v as 
opposed to the propagation rate of a particle. Substitution from equa- 
tions (36) and (37) reduces this to 


Dy 

Dt 





du 

dt 


or with the aid of equations (2.6), 


(27), (29), and (30) 


Ik 



Note that Dy/Dt may be interpreted as the "wave velocity” of the portion of 
a disturbance with velocity v. Denoting this wave velocity by c and substi- 
tuting from equations (20) and (3l)> one obtains 


c(A) 




(40) 


Not only does v propagate with velocity c(A) but so does A; 3 therefore 
equation (40) is the relation required to calculate the distortion of a large - 
amplitude wave as it propagates into the undisturbed tube. The only quantities 
required are (l) the pressure-area curve for the tube f(A) , (2) the density 
of the fluid, p, (3) the undisturbed cross-sectional area of the tube, A e , and 
(4) the spatial distribution of the cross-sectional area (or pressure since 
there is a one-to-one correspondence by equation (12)) at some initial time 
t Q , A(x,t 0 )« The actual calculation of the distortion will be considered in a 
later section. 


Criteria for Steepening 

Of primary concern in considerations which follow will be the waves which 
steepen as they propagate. A steepening wave is defined as one in which the 
absolute value of the slope (|Sp/Sx| or |Sa/Sx|) for some value of pressure. 


3 The derivation of the propagation velocity of A differs from that of 
v only in the determination of the particle -to -particle propagation rate. 
Instead of equation (35) > one obtains 


dA(x,t) + SA(x,t) 
Sx St 


and with equations (20), (2l), and ( 27 ) 


Sa 

Sx St 

st Sa 

Sx 


Se(x,t ) 
St 

Se(x,t) 

Sx 



which is identical to equation (36). Since the pressure p is assumed to 
have a one-to-one correspondence with the area A, then the function c(p) can 
be calculated from c(A) . 



or cross-sectional area, increases as the -wave propagates- If a steepening 
wave propagates for a sufficiently long time a sharp wave (infinite slope) 
will develop (at least mathematically). 


One could consider the steepening phenomenon by evaluating dc/dA, that 
is, from equation (40) 


dc(A) 

dA 


_ 1 
+ o 


f y (a)a' 


“ 1/2 


p[f 1 (A) 4- Af "(A) ] + 


" f f (Ay 

Ap 


, 1 / 2 


m 

\ 


However, equation (4l) does not yield a simple set of criteria for steepening 
waves due to the complexity of the right-hand side of the equation. However, 
if c is written as a function of the parameter a (eq. ( 32 )) one obtains 
with the aid of equations (33) and (39) 


c(a) 



2a) f "(a) 
P 


+ 



f y (a) 


p(l - 2a) 


1/2 


da 


(42) 


Differentiation of 


c(a) in equation (42) produces 


dc(a) 

da 



.1/2 


1 - 2a 
pf *(a)_ 


f"(a) 


(43) 


Recall that the positive sign in equation (43) is applicable to -waves 
traveling to the right (positive x) and the negative sign applicable to waves 
traveling to the left (negative x) . Therefore, for waves traveling to the 
right 


dc 

da 

0 

if 

f "(a) 

dc 

da 

0 

if 

f"(a) 

dc . 
da 

0 

if 

f"(a) 


Rote also that for compressive waves 


(44) 


A e „ 

— decreases, a increases 


an d for decompression waves 


A e . 

-7- increases , 


a decreases 


(45) 

(46) 


l6 



The criteria for steepening -waves may now Toe stated in terms of the 
pressure as a function of the parameter a. The most useful criteria are 
probably the ones in terms of a plot p(a). 

If one plots the curve p[l - (A e /A) 2 ] or p(a), the criteria for steepen- 
ing waves may be stated as: 

(a) Compression waves will steepen for ranges of p in which the slope 
of the p(a) curve is increasing (f"(a) > 0) . Conversely, they will flatten 
otherwise. 

(b) Decompression waves will steepen for ranges of p in which the slope 
of the p(a) curve is decreasing (f"(cc) < 0) . Conversely, they will flatten 
otherwise. 

(c) Compression and decompression waves will propagate undistorted for 
ranges of p in which the slope of the p(a) curve is constant (f M (a) = 0) . 

The last of the criteria is in accordance with the condition of nondis- 
torting waves (eq. (19))- It seems appropriate to reemphasize that the 
criteria for steepening are stated entirely in terms of a function relating 
the. pressure to the area of the tube. 


II. CALCULATION OF WAVE DISTORTION 


General Equation of Distortion 

The purpose of this section is to establish the formulas for calculating 
the distortion of a large -amplitude wave propagating into an undisturbed tube. 
In addition, later parts of this section will be devoted to determining the 
criteria for a sharp wave (infinite slope) and the critical distance required 
for a sharp wave to develop. Finally, an example calculation will be presented 
to demonstrate an application of the theory. 

Section I provides the theory necessary to determine the spatial wave 
velocity as a function of the pressure (or area) amplitude. If this quantity 
is denoted by c(p), then it follows from the discussion of the first section 
that the pressure, p, maybe expressed as a function of t - [x/c(p)], namely 


p(x,t) = h 


t 


X 

c(p) 


( 47 ) 


where c(p) is a known function. Equation (47) is applicable to a wave 
traveling to the right (positive x) . 
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Sharp-Wave Criterion 


The pressure gradient will become infinite when 

- oo (48) 

dx 

Evaluating dp/dx from equation (47) yields 


dp(x,t) _ 

Bx 


x 

c(p)„ 


-c(p) 


+ X 


Mv) 

<^x 


[c(p )] 2 


but 


dc(p) _ dc(p) Sp 
8x <3.p 8x 


therefore > 


dp(x,t) 

Sx 


t 


X 

c(p). 


1 X dc(p) dp(x,t) 

c(p) C 2 (p) <3-P Sx 


or, solving for dp/dx , 


dp(x,t) 

Sx 


h f 


t 


x 

c(p). 


c(p) 


dc(p) 1 
dp c(p) 


b 1 


t 


X 

c(p) 


( 49 ) 


The function b[t - x/c(p)] is assumed to be known for some tube position, 
namely, x = 0, so that h(t) is a prescribed function. Since this study is 
concerned with the transition from "smooth" waves to sharp waves, h(t) is 
assumed to be continuous and dh(t)/dt (or h*(t)) is assumed to be finite. 

The pressure gradient 3p(x,t)/8x becomes infinite if 


1 


x 


dc(p) 

dp 


c 2 (p) 


h* 


t 



(50) 


The criterion for a sharp wave to develop (eq. (50)) depends on (a) the pre- 
scribed form of the pressure, h(t), (b) the wave velocity pressure curve, 
c(p), and (c) the location on the tube, x. 
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Critical Length 


the 

if 


As the pressure wave passes some position x = L c along the tube, 
slope of the wave Sp(L,t)/Sx will he infinite for pressure level p* 
(from eq. ( 50 )) 


L c = 


e 2 (p*) 

dc(p*)/dp 



(51) 


Note that t Q is the time at which p(Lc,t) = p* which is (from eq. (47)) 


p*(L c ,t 0 ) = h t 0 - 


Lc 


:(P*). 


or 


to = -T~T + L _1 Ep*] (52) 
c(p*) 

where h _1 denotes the inverse function and should not he confused with h*, 
which denotes a derivative with respect to the argument of the function, nor 
with the negative one power. The denominator of equation (5l) can he trans- 
formed hy substitution from equation ( 52 ) to obtain 


h» 




h' [h -1 (p*) ] 


(53) 


Differentiation of both sides of equation (47) with respect to t yields 


h' 


t 


x 

c(p)_ 


dp(x,t) 

St 

x[dc(p)/dp] [Sp(x,t)/St] 
c(p) 2 


(54) 


The right side of equation ( 53 ) > h f [h -:L (p*) ] , is the value of h 1 evaluated 
for the value of the argument at which h = p*. This may he evaluated from 
equation (5*0 "by choosing x = 0 and evaluating the right side of the equation 
for p = p* thus 


h*[h 1 (p*)] 


h*(t) 


h=p* 


Sp(0,t) 

St 


p=p* 


Sh(t) 

St 


h-p* 


(55) 


where [Sh(t)/St] 1^-^* can he readily calculated from the known function 
h(t). P 
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the compression waves sharpen -while 
the decompression waves become less 
steep so that at some time ti the 
wave has the shape shown in fig- 
ure 5 * The pressure in the tube is 


y 


p(x,t) = Pi sin (A) 


t 


c(p)J 

(60) 


The computation of the critical 

Figure 5.- Spatial distortion of pressure wave length Lc(p*) requires the evalua- 

for illustrative example. tionof dh/dt | h=p *. Krom equa- 

tion (58) , 


= Pild COS lot 

ot 


and the time at which h = p* is 


t 


h=p 


* 


1 . -1 

u sin 


Pi 


Therefore 


dh 

St 


= PilO cos 

r . -1 p*i 

sin — 

h=p* 

P!_ 


or simply 


= ± (0 V Pi - p* 2 ( 61) 

h=p* 

where the choice of ± sign depends on the position on the function h(t) 
and, therefore, the position on the wave. 

The numerator term in the calculation of L c (p*) is, from equations (56) 
and ( 59 ), 


6h 

3 t 


c* 2 ^*) 

dc(p*)/dp 


cipi exp 



(62) 


Substitution from equations (6l) and (62) into the critical length 
equation yields 


Lc(p *) = cip! exp(pVpi ) 

± 10 n/ p 2 - ' p * 2 


21 



or 


Mp*) 


(63) 


Cl exp(p*/p 3 .) 

±w \/l - (p*/pi) 2 


In order to obtain the minimum critical distance,, differentiate the right side 
of equation (63) with respect to p* and equate the result to zero 


aL c (p*) 

dp* 


ci exp(p*/p!) 1 

pi w ±[1 - (p*/pi) 2 i 3/2 L pi 



(6k) 


thus a relative minimum or maximum of the function 


L c (p*) occurs when 


pf 1 -{5. 

Pi 2 + 2 


but the range of interest of p*/pi is -1 5 p*/pi £ 1. Therefore, the 
minimum of interest is when 


Pi 


2 “ -0.6l8 


or from equation (63) 


L cm = — °- 686 


The negative pressure is a differential pressure in the previous relations. 
The absolute pressure must, of course, be positive. The computation of the 
minimum length of tube required for the development of a steep wave is 
extremely simple for the example problem. Although this simplicity will not 
prevail for the more general problem (i.e., more complex c(p) and h(t) 
functions ) , a numerical calculation should always be practical. 


III. SHOCK WAVES 


The formulas, developed for calculating the length of tube required for a 
sharp wave (dp/6x = 00 ), strictly speaking, do not apply in the limit 
(dp/chc = 00 ) but should be applied only as 8p/8x approaches infinity, for 
when dp/dx = 00 the particle acceleration also becomes infinite. Also, for 
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very sharp pressure gradients the axial bending rigidity of the tube wall 
should be included in the analysis. 4 

The final transition from the nearly sharp to a truly sharp wave presents 
a formidable mathematical problem. However, the propagation of a sharp wave, 
once the transition phase has been completed, can be treated quite simply 
provided one again assumes that the axial bending rigidity of the tube is 
negligible. 

Associated with the sharp waves are discontinuities in the particle veloc 
ity and pressure and (except for the linear wave equation) energy losses. For 
the remainder of this report the sharp waves will be called "shock waves." 


Shock -Wave Velocity 


Consider the propagation of a compression shock wave through a fluid- 
filled elastic tube (fig. 6). Let the cross-section area, the pressure, and 

the particle velocity (A, p, and v, 
respectively) be distinguished by sub- 
scripts 1 and 2 referring to the right 
and left side of the shock wave, respec 
- tively. The spatial velocity of the 
shock wave is denoted by eg . Conserva 
tion of fluid mass across the shock 
front requires that 


^2 * p 2 ’ 


A h P|» v l 


Figure 6.- Shock-wave notation. 


Vi)Ai = (c s - v 2 )A 2 ( 65 ) 


A useful parameter in further computations is c* which is defined as 
the shock-wave velocity in the compressed (or lengthened) coordinate system. 
The relation of cs to Cg may be written 


c s - vi + c|*( €i + l) 


( 66 ) 


where e x = A e /Ax - 1 (see eq. (20)). Or if one had chosen to write the rela- 
tion in terms of the left side of the shock 


c g - v 2 + c*( e 2 + l) 


(67) 


4 The axial bending rigidity in the development of a shock requires fur- 
ther investigation. This effect may or may not be important. Certainly, how- 
ever, it would be more attractive from the physical viewpoint if the tube wall 
displacement discontinuity were eliminated. Indeed it is not certain that a 
shock wave will develop if the bending rigidity is included. 
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Note that c g of equations (66) and ( 6 7) must he equal, therefore 


(v 2 - Vi) - c*( e i - € 2 ) 


(68) 


which is equivalent to equation ( 65 ). 

Besides mass, the momentum of the fluid must also he conserved. The 
impulse of the forces on the fluid acting along the tube is 


P2A2 - P1A1 


oA 2 

Ai 


f(A)dA 


Ax 

c s “ vi 


which must equal the change of momentum of the fluid 


pAxAi(v2 - vi) 


therefore 


I-P2A2 - P1A1 

Elimination of c s and V 2 
and (68) produces 



f(A)dA 



from equation ( 69 ) with 


= pAx(v 2 - vi) 
the aid of equations 


(69) 

(66) 


or 


( c *) 2 = — — 

V t 'S' ^A. 


pAi TTTTJ (T x - 62) 


P2A2 “ P1A1 



f(A)dA 



(70) 


Integration by parts of the numerator of the right-hand side of equation (70) 
simplifies the expression for c* 


( 4) 2 = 

and the shock wave velocity in the 
eqs . (66) and (7l)) 


x A 2 


An 


A Mihl 

dA aA 


Ae A e 


pAe \Ai " A 2/ 
undeformed coordinate system is (from 


(71) 
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a 


1 

pAf 


r A2 A u 
J Ai dA 


( 72 ) 


(c 


s 


- V3.) 2 = 


(4Y 




The velocity of the shock wave can be determined if the pres sure -area 
curve of the tube, the conditions ahead of the shock ( A x , vi), and the 
amplitude of the shock A 2 - Ai (or p 2 - Pi) are prescribed. 


The next step is to estimate the shock -wave velocity in terms of the 
small -amplitude wave velocities of disturbances at the foot (A = Ai) or peak 
(A = A 2 ) of the original wave (fig. 7 ) • This is accomplished by introducing 




Smooth wove Shock wave 


Figure 7.- The pressure as a function of a for use in shock -wave description. 
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1 


the a notation (eq. ( 32 )) into equation (7l) 

P a 2 df(a) 


(c*) 2 


a-i 


da 


7 1 - 2a 


da 


) (7 1 - 2ai - 7l - 2a 2 ) 


(73) 


or from equation ( 72 ) 


(1 


(c s - vi) 2 = 


oa 2 

2a x ) / 

^ ai 


df(a) 


da 


7l - 2a 


da 


)(7"l - 2ai - 7 1 - 2a 2 ) 


(74) 


The equations derived thus far in this section are valid for infinitesimal as 
well as finite wave amplitudes. The wave velocity of the foot of the wave 
which developed into the shock wave was (from eq. (74)), taking the limit as 
a-2 -* ax, 


( c s - 

and the velocity at the peak of 

( c s " 


Vl)f = 


1 - 2ai df(ai) 


da 


the original wave^ was 
.2 1 - 2a 2 df(a 2 ) 


vi 


)n - 


da 


(75) 


(76) 


Note that if df/da is positive and d 2 f/da e is positive for a x ^ a £ a 2 
(i.e., df(a x )/da < df(a 2 )/da) as in figure 7 (which is the criterion for sharp 
waves, section i) , then from equations (74), (75) and ( 76 ) 


(c S " Vl)f < (c s - vi) < (c s - Vl)i 


(77) 


The velocity of the developed shock wave is greater than the velocity of the 
foot of the original wave hut less than the velocity of the peak of the 
original wave. In the limiting case in which df/da = constant (nondistorting 
wave, see section i) then 


( c s _ 


vi)f = 


(c s - Vx) 2 = (c s 


- Vi)§ 


(78) 


5 Strictly speaking, the "original" wave is a smooth wave with amplitude 
equal to the amplitude of the shock wave. It has not heen proved that the 
smooth wave will develop a shock wave of equal amplitude. 
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Energy Loss 


Since the shock -wave velocity is determinable from considerations of 
conservation of mass and momentum, the energy loss across the shock can be 
calculated directly. 

The external work rate on the fluid is 

(T9) 


( 80 ) 

the rate of change in potential energy of the tube is 


P2A2V2 - P1A1V1 

the rate of change of kinetic energy of the fluid is 

\ pAi(v| - vf)(c s - vi) 


■' A 2 


f (A)dA Cg 


and the energy loss rate is defined as 


( 81 ) 


(c s - Vi)El 


( 82 ) 


Equating the work rate to the energy change rate yields 

1 P A 2 

P 2 A 2 V 2 - P 1 A 1 V 1 = p pA 2 (vI - vf)(c s - vi) + / f(A)dAc s + E L (c s - vi) 

J Ai 

(83) 

If V 2 and c s are eliminated from equation ( 83 ) with the aid of equations ( 65 ) 
and ( 72 ) then, after considerable manipulation of terms, one obtains 




p 2 (A 2 - Ai) + PiA x 


-it 1 - 


1 + 


it) L 


a 2 


f(A)cLA 


(84) 


Thus the energy loss coefficient Ep can be calculated if the conditions 
ahead of the shock wave, the amplitude of the shock wave A 2 - Ai (or 
p 2 “ Pi)j> and the pressure-area curve for the tube are known. 

The case considered here is one of a particular class of shocks which 
require only the mechanical shock conditions (ref. 10) for the determination 
of the shock. The energy loss is generally attributed to the generation of an 
equivalent amount of heat. 
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IV. KOROTKOFF SOUNDS 


Sphygmoman ome try 

Although blood pressure can be measured directly by inserting pressure 
measuring devices into the arteries or veins , indirect methods of measurement 
have been devised which provide less discomfort to the patient, eliminate the 
risk of infection, and provide greater convenience to the medical practitioner. 

Several indirect methods (ref. 11 ) have been devised and used success- 
fully; however, one method, sphygmomanometry, has come to be accepted as a 
standard in the medical field. Fortunately, the method is simpler than its 
name and is familiar to almost everyone as part of a physical examination. 

The details of the method may not be so familiar to those outside the medical 
profession and, at the risk of being superfluous, will be presented. 

The required equipment consists of a compression cuff, an inflation bulb 
with pressure control valve, and a mercury manometer or other pressure measur- 
ing device. 'When the cuff is snugly applied to the arm just above the elbow 
(fig. 8), inflation of the cuff with the bulb compresses the tissues which in 



turn compress the brachial artery. The pressure in the cuff is increased 
until no pulse wave can be felt in the artery downstream from the cuff. 

A stethoscope is positioned just downstream from the pressure cuff and 
the pressure in the cuff is slowly reduced by the pressure control valve. As 
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the cuff pressure decreases, certain distinctive sounds may he heard until the 
cuff pressure is reduced to a sufficiently low level at which the sounds disap- 
pear. The sounds are called Korotkoff sounds and the pressures indicated on 
the manometer (cuff pressures) at which the sounds appear and disappear are 
generally accepted as a measure of maximum blood pressure (systolic) and 
minimum blood pressure (diastolic) in the brachial artery. 


Direct and Indirect Method Correlation 

The simplest and most common correlation of the sounds and the actual 
pressure pulse is shown graphically in figure 8. This correlation assumes 
that the highest pressure at which sounds first occur corresponds to the maxi- 
m-urn value of the pulse pressure and that the lowest pressure at which sounds 
disappear corresponds to the minimum value of the actual pulse pressure. 

Several medical researchers have attempted to establish the validity of 
the indirect method by comparison with blood pressure measurements obtained 
from intra-arterial pressure transducers (refs. 12-15). Although the degree 
of correlation varies between investigations (the medical history of the sub- 
jects varies also), it is not unusual to obtain discrepancies between the two 
methods as great as 20 percent. There are, of course, cases when the indirect 
method produces zero diastolic readings - a grave condition if the method were 
correct. 

It seems unreasonable to expect that the discrepancies between the two 
methods will be understood until the conditions required for the generation of 
the Korotkoff sounds and the source of the sounds are known. 


Various Hypotheses 

Hypotheses explaining the Korotkoff sounds have been plentiful. 6 No 
attempt will be made to be all inclusive but four of the more popular 
hypotheses from the literature will be briefly reviewed. 

Water hammer . - Early in his investigation of blood pressure estimation by 
indirect methods, Erlanger (ref. 17) (1916) came to the conclusion that the 
main mechanism of Korotkoff sounds is as follows: 

"Under compressions which permit the pulse to determine relatively 
wide excursions of the arterial wall in the compression chamber, 
that is, under compressing pressures ranging from systolic arterial 
pressure to, and even a variable distance below, diastolic pressure, 
the volume of the compressed artery increases abruptly with each 
pulse. This permits a considerable volume of blood to enter the 
opening artery with a high velocity. The motion of this column of 
blood is, however, suddenly checked where it comes into contact 
with the stationary, or practically stationary, column of blood 
filling the uncompressed artery below. The water hammer that is 

6 The introductions to references 1 6 and 8 provide historical reviews. 
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thus set into play distends the arterial vail at the point of impact 
with unusual violence. This distention sets the arterial vail into 
vibration and the sound is produced.” 

After further experiments, however, Erlanger (ref. 18) changed his 
opinion and supported the hypothesis of Bramvell. 

"Breaker 11 phenomenon . - While vorking in Erlanger f s laboratory, Bramvell 
(refs. 19 and 20 ) ( 1923 ) experimentally studied the change in the form of the 
pulse vave in the course of transmission. He correlated the propagation of 
the pulse vave vith the propagation of ocean vaves tovard a beach and 
concluded that for the pulse vave: 

"Hence, the more rapidly moving elements vill tend to catch up 
to those vhich are moving more slowly; in other vords the vave vill 
tend to become more vertical, and, under suitable conditions, parts 
of it may actually tend to topple over and form breakers . * Even if 
the vave as a vhole does not *break, f any individual portion of it 
vhich does so vill set up irregular vibrations in the arteries, 
vhich, if they be of sufficient amplitude, vill be transmitted to 
the periphery . . . . " 

Later Bramvell related the "breaker" to the sound production in arteries 
and specifically to the Korotkoff sounds. 

Bernoulli effect . - Rodbard (ref. 21 ) experimentally studied the flov of 
fluid through collapsible tubes. Due to the coupling of the flov (through the 
Bernoulli or flov velocity pressure effect) and the flexibility of the tube, 
he vas able to produce instabilities and vibrations of the fluid -filled tube. 
He termed this phenomenon "flitter" and one of his conclusions vas 

"The Korotkoff murmurs heard vhen arteries are compressed during the 
measurement of blood pressure are probably also due to flitter of 
the vails . . . . " 

It should be noted that the "flitter" of vhich Rodbard speaks requires 
relatively high flov velocities and it is known that the Korotkoff sounds can 
be reproduced vith zero mean flov (ref. 22). 

Dynamic instability . - Anliker and Raman formulated a hypothesis for the 
Korotkoff sounds at diastole (ref. 1 6 ) and Raman later extended the hypothesis 
to include systole (ref. 23 ). The hypotheses interpret the sounds as "a phe- 
nomenon of dynamic instability (oscillations vith increasing amplitude), the 
instability being induced by the application of a pressure cuff." They 
explain their hypothesis as follovs: 

"Like any other sound perceived by a human ear, the Korotkoff 
sounds heard vith the aid of the stethoscope are aggregates of vibra- 
tions vhose frequencies and intensity place them vithin the audibil- 
ity range. It is conceivable that the Korotkoff sounds are due to 
disturbances in the flov that are induced locally, that is, in the 
segment of the brachial artery vhich is compressed by the cuff. 
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However , it is also possible that they are a part of* the noise asso- 
ciated with pulsative flow through the circulatory system that is 
selectively amplified to a level above the audibility threshold in 
the compressed section of the brachial. It is well known that pul- 
sating flow through the complex system of blood vessels generates a 
multitude of wave motions; the associated vibrations are not audible 
because their amplitudes and in many cases their frequencies are far 
below the audibility threshold. Irrespective of the origin of the 
disturbances leading to Korotkoff sounds, it will be shown that the 
application of a pressure cuff can change the system locally from 
one that is dynamically stable to one that is intermittently unsta- 
ble and thus capable of mechanically amplifying certain low 
intensity vibrations (disturbances) to such an extent that they 
become audible . . . . M 

Theoretical work to support the hypotheses is practically nonexistent 
except for that of Anliker and Raman. They have done the linear analysis for 
a fluid-filled circular cylinder (simulating diastole) and Raman did a similar 
analysis for systole except that the shell geometry was changed to simulate 
the "flattened” artery. Their analyses do not provide a means of computing 
the frequencies or amplitudes of the sounds (which are assumed to be amplifica- 
tions of disturbances already present) but establishes the differential 
pressure for instability of the system (buckling of shell) and amplification 
of the disturbances. 

Bramwell and Hill (ref. 19 ) did make some calculations of the distortion 
of the arterial pulse wave assuming that the small displacement wave velocity 
of the pulse wave was known. Their calculations were, however, incorrect 
because they, in effect, included only one term of the right side of equa- 
tion (40). In particular, they used only the small displacement velocity 
(measured) which corresponds to 


f f (A)A 


1/ 2 


P 


They have therefore neglected the term 



1/2 


dA 


or the "integrated particle velocity." 


Applicability of Present Theory 

The intent of this section is to investigate the fundamental physical 
phenomena associated with the production of Korotkoff sounds. Parameters that 
influence the character of the sounds, but are nonessential for production of 
the sounds, will be neglected. This is done to isolate the physical features 
essential to the production of the sounds. Future analysis should provide a 
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method for determining the effect of those parameters which influence but are 
not essential to the production of the sounds. 

The applicability of the analysis of sections I through III to the inves- 
tigation of the Korotkoff sounds depends on the ability of the mathematical 
model to represent the important properties of the physical system which 
produces the sounds. This analysis will be applicable if it can be 
demonstrated that: 

(1) The phenomenon of Korotkoff sounds can be simulated by pulsatile flow 
in a compressed isolated artery (that is , an artery detached from its 
surrounding tissue but still intact at its ends), and 

(2) The assumptions of section I are applicable in formulating the 
mathematical model for the pulsatile flow. 

The validity of the first condition has been demonstrated in experiments 
on anesthetized dogs by Erlanger (ref. 2k) and in experiments on simulated 
arteries by Sacks et al. (ref. 25). The results of these experiments demon- 
strate that the Korotkoff sounds occur without the surrounding tissue. This 
should not be taken to imply that the surrounding tissue will not change the 
quality of the sounds. Anliker and Raman (ref. 1 6), in their experiments on a 
simulated artery, have shown that simulated tissue surrounding the artery 
muffled the sounds. The quality of the sounds is affected by the tissue but 
the tissue is not essential to produce the sounds. 

The second condition to applying the analysis will be satisfied by 
considering each assumption of section I explicitly. 

(a) The effect of compressibility of the blood in the study of pulsatile 
flow in the arteries has been shown to be small (ref. l) . The effect of vis- 
cosity on the Korotkoff sound phenomenon was considered experimentally by 
Sacks et al. (ref. 26). They concluded that the effects of viscosity are rela- 
tively small since it was demonstrated that Korotkoff sounds could be produced 
over a range of viscosities from 1 to 30 cm 2 /sec. The kinematic viscosity of 
blood is of the order 5 cm 2 /sec. 

(b) The assumption of one-dimensional flow and the neglection of radial 
inertia have been shown to be reasonable approximations as a result of 
experiments performed by the author to supplement the available experimental 
data. These experiments are discussed below. 

In the author’s experiments the simulated artery was either a gum rubber 
tube ( 1-inch diameter and l/l6-inch wall thickness) or a transparent flexible 
plastic tube of similar dimensions. The fluid was either pure water or a 
water-glycerin mixture (90-percent water, 10-percent glycerin). The pulsatile 
flow was provided by a positive displacement pump (fig. 9)* The mean flow was 
zero in accordance with the findings of Sacks et al. (ref. 25) that the mean 
flow effect was negligible. The compression (simulated cuff) of the tube was 
provided by decreasing the mean internal pressure and thus allowing the 
atmospheric pressure to provide the cuff pressure. This method of compressing 


32 



- 2 ' 


Rigid tube reducer and 
transition zone 



f 


I diameter 


2" diam piston driven 
sinusoidally with Scotch 
yoke 


z 



Reservoir system (closed) 
50' f f- 


Stethoscope 

^Displacement 
\transducer (typical) 



Flexible tubing supported 
on rigid base 


Figure 9*- Diagram of test fixture for author's experiments. 

the tube eliminates any influence a mechanical cuff may have on the sounds. 7 
In addition, this method allows almost unlimited access to the tube for 
instrumentation. 


Flow visualization experiments were performed using the clear plastic 
tube and a solution of water (9° percent) and glycerin (lO percent) of the 
appropriate density to suspend small plastic beads (approximately 0 . 05 -inch 
diameter). With the pump providing a sinusoidal displacement input at the end 
of the tube the mean pressure differential across the tube wall was adjusted 
to produce sounds at the stethoscope location. High-speed motion pictures 
were then taken of the area surrounding the stethoscope so that the tube and 
fluid motion could be studied in detail. The pictures showed that the veloc- 
ity of the suspended particles and thus the velocity of the fluid was uniform 
across the cross section of the tube. In addition, the maximum lateral motion 
of the particles was about 7 percent of the longitudinal motion. On the basis 
of these observations the assumptions of a one -dimensional flow and negligible 
energy associated with the lateral motion are justified. 

(c) Since the density of the tube wall is approximately the same as the 
density of the fluid, the mass of the tube can be neglected for the same 
reason the lateral inertia of the fluid can be neglected. 

(d) Although prestress of the tube in the axial direction affects the 
quality of the Korotkoff sounds (ref. l8), prestress is not essential to the 
production of sounds. This was demonstrated experimentally by Bramwell and 
confirmed by the author in his experiments. 

The flexural rigidity in the axial direction of the tube has no effect 
for long wave lengths (ref. 3) but the effect must be accounted for as the 
pressure gradient becomes sharp. The transition region (pressure gradient) at 
which the flexural rigidity must be included requires further investigation. 
The justification for omitting axial flexural rigidity in the present 

T The "cuff" in the author 1 s experiments was actually the same length as 
the elastic tube (approximately 50 feet). The "effective" cuff length was, 
however, determined by the location of the stethoscope. If the pulse wave 
form were measured at station xi and the stethoscope were located at station 
X 2 (downstream from xi) then the effective length of the cuff was x^ - xi. 
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investigation of Korotkoff sounds is based on the good quantitative agreement 
of the author f s analysis and the experimental data of Bramwell (presented in a 
later paragraph) . 

(e) The assumption of uniformity of the tube under equilibrium conditions 
is certainly justified for the simulated artery experiments. The application 
to the brachial artery is validated by the work of Tickner and Sacks who made 
extensive measurements on excised brachial and other arteries (ref. 27 ). 


Hypothesis for Mechanism 

The present hypothesis for the Korotkoff sounds is divided into two parts. 
The first considers the development of necessary conditions (mechanism) for 
the production of sounds and the second part considers the sounds. The 
proposed mechanism should already be apparent from the development of 
sections I through III. The hypothesis is: 

As the pulse wave traverses the compressed portion of the 
artery, it may, with favorable system parameters, steepen and, in 
fact, become sharp (possibly a shock wave). The system parameters 
include (l) the shape of the arterial pulse as it enters the com- 
pressed section of the cuff, (2) the physical properties of the tube 
wall, (3) the cuff pressure, and (4) the length of the cuff. The 
development of a shock wave is not essential to the production of 
sounds but the wave must become relatively steep (see the sound 
hypothesis) . 

It should be noted that the cuff pressure could be zero and the other 
parameters favorable to produce steep waves. This could account for the zero 
diastolic reading sometimes obtained in practice. 

This hypothesis does not differ significantly from that of Bramwell 
except for the final stage of the steepening wave which is assumed here to be 
a shock wave rather than the less probable "breaking wave." The major contri- 
butions to this hypothesis are (i) the theoretical support of sections I to 
III which make the hypothesis both plausible and useful, and (ii) the 
explanation for the developed sharp wave. 


The Sounds 

The mechanism hypothesis offers no explanation for the sounds heard 
through the stethoscope. Generally, one associates sound with oscillatory 
pressures. However, no oscillatory components of pressure are predicted by 
the analysis and in the author 1 s experiments sounds could be obtained when the 
displacement of the tube wall contained no oscillatory components. High-speed 
motion pictures (2000 frames per second) and tube wall displacement trans- 
ducers revealed only the sharply rising and slowly decaying wall displacement 
as each "pulse" wave passed. In addition, the flow visualization studies 
showed no turbulence in the flow when the sounds were heard or during any 
other part of the cycle. This last observation is in agreement with the 



turbulence studies of Anliker 
and Raman (ref. l6) who used 
dye for flow visualization. 
Measurements from pressure 
transducers inside the tube 
(fig. 10) showed no oscilla- 
tory pressure components , only 
the rapid pressure rise as the 
"pulse” wave passed the 
transducer. 

The author f s experimental 
results and the analysis indi- 
cate that the predominant 
source of the sounds is the 
sharp pressure gradient. In 
order to determine if this 
alone were a plausible explana- 
tion^ a series of additional 
experiments were conducted. 

The signal from a pres- 
sure transducer inside the 
tube and flush with the inner 
wall of the tube was used as the "driving" signal for an audio amplifier 
speaker system (fig. ll) . A stethoscope was placed on the tube wall directly 
over the pressure transducer. The sounds produced by the speaker system were 
then compared with those obtained by listening through the stethoscope. The 
comparison of sounds by the human ear is at best qualitative; however, the 
sounds from the speaker and those from the stethoscope were essentially but 
not indistinguishably the same. In order to insure that the speaker sounds 
were not affected by artifacts of the pulse pressure wave which were not 
observed on an oscilloscope display, a function generator was used to simulate 
the pressure transducer waveform. The function chosen was a ramp with the 
slope approximating the pulse pressure gradient. When this signal from the 
function generator was Introduced into the speaker system the sounds were 
indistinguishable from those obtained when the pulse pressure signal was used. 
It was concluded that the sharp pulse wave alone could account for the 
predominant sounds. 



Figure 10. - Typical pulse wave and simulated Korotkoff 
sounds. Upper trace is pressure recorded from 
pressure transducer inside tube under stethoscope 
(see fig. ll) . Lower trace is signal from microphone 
mounted inside stethoscope earpiece. Time base 
0.2 sec per division. 



Figure 11.- Relative location of stethoscope and pressure transducer. 

As noted above, the sounds from the stethoscope were "essentially" but 
not "Indistinguishably" the same as those from the speaker. Since a 
stethoscope has resonant frequencies ("organ pipe" frequencies) that result 
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from the tubing connecting the head of the stethoscope to the earpiece, it was 
decided to investigate the possibility that the "difference" of sounds could 
be due to the stethoscope resonances. This study was made by varying the 
length of tubing connecting the stethoscope head- and earpiece. The variation 
of sounds as a function of tube length was observed by two methods - listening 
through the earpiece and using a microphone inside the earpiece so that the 
sound pressure waveforms could be observed on an oscilloscope (fig. 12). The 
oscilloscope traces showed the expected trend of increasing frequency in 
pressure fluctuations as the stethoscope tube length was decreased. The 
sounds appeared to become correspondingly "sharper" as the length of 



(a) 11 feet of tubing. 



(b) 5-5 feet of tubing. 


Figure 12.- Typical pulse wave and simulated Korotkoff sounds for various lengths of tubing connecting 
stethoscope head- and earpiece. Upper trace is pressure recorded from pressure transducer inside 
tube under stethoscope (see fig. ll) . Lower trace is signal from microphone mounted inside 
stethoscope tubing. Time base 0.1 sec per division. 
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(c) 2 feet (standard) of tubing. 



(d) 0.4 feet of tubing. 
Figure 12.- Concluded. 


stethoscope tubing was decreased and more closely approximated the speaker 
system sounds (from the pressure pulse wave). The author concluded that the 
differences in the sounds from the stethoscope and the audio system (with 
pulse wave as input) were accounted for by the characteristics (resonances) of 
the stethoscope. 
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The range of frequencies 8 of the Korotkoff sounds reported in the litera- 
ture is quite varied and includes the threshold of hearing (approximately 
20 cps) to about 500 cps (refs. 8 , 22, 30, and 3l) • During the development of 
an automatic hlood pressure monitoring system for the Mercury manned space 
flights (ref. 29 ) it -was discovered that the microphone signal containing the 
Korotkoff sounds could be filtered with negligible loss of accuracy in blood 
pressure determination. It was determined that if the frequency band pass 
were set so that only components between 32 and 40 cps were passed, the indi- 
rect blood pressure determination was essentially the same as if the unfil- 
tered signal were used. This result implies either that the predominant 
oscillatory components of the sounds were in the frequency band of the filter 
(32 to 40 cps) or that the microphone filter was not "filtering" the sounds 
but was "responding" to a sharp input. In light of the previous discussion 
the latter explanation seemed more plausible to the author and a test was 
conducted to verify this suspicion. 

The pressure transducer signal was filtered by a band-pass filter and 
then introduced into the audio speaker system. The pressure waveform before 
and after filtering was monitored on an oscilloscope (fig. 13)* It was found 
that the filter changed the input waveform from a sharp rise and slow decay to 
a momentary decayed oscillation. The sounds were not eliminated but their 
quality changed so that they did not seem as "sharp" as the unfiltered sounds. 
The frequency of the filtered signal was that of the filter (about 3 8 cps when 
the bandwidth was 3 2 to 40 cps). In order to confirm that the oscillatory 
decay was due to the response of the filter (and not the frequency components 



(a) 12-16 cps "band-pass filter. 


Figure 13.- Typical pulse wave and pulse wave after narrow "band filtering. Upper trace is 
pressure recorded from pressure transducer inside tube. Lower trace is same signal after 
narrow "band filtering. Time base 0.2 sec per division. 

8 The frequencies referred to are those of the oscillatory components of 
the Korotkoff sounds when displayed on an oscilloscope or other recording 
device (ref. 28, for example) and should not be confused with the Fourier 
components of a pulse wave type record as reported by Ware (ref. 29 ). 
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(b) 22-26 cps band-pass filter. 



(c) 32-36 cps band-pass filter. 
Figure 13 *' Concluded. 
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of the Korotkoff sounds) a square wave function generator was used in place of 
the pulse pressure. The filtered response once again was a series of short 
hursts of decayed oscillations with frequencies of the filter "center” 
frequency (fig. 14) . 



(a) 12-16 cps band-pass filter. 



(b) 22-26 cps band-pass filter. 


Figure 1^.- Square wave and narrow band- pass filtered square wave. Time base 0.2 sec per 

division. Upper trace is square wave from function generator. Lower trace is same signal 
after narrow band-pass filtering. 
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(c) 32-36 cps band-pass filter. 

Figure 14. - Concluded. 

Much work remains to he done to confirm that all the sounds observed in 
sphygmomanometry are due to the phenomenon of the sharp wave. The author 
believes that the preliminary experiments and analysis presented in this 
report offer sufficient evidence for proposing a hypothesis for the sounds. 
This is, of course, only a hypothesis that must be evaluated by further 
investigation both analytical and experimental. 


Hypothesis for Sounds 

On the basis of the analysis for the mechanism of the Korotkoff sounds 
presented previously and experiments by the author, it is concluded that: 

The Korotkoff sounds are the response of the measuring system to 
a steep (possibly sharp or shock) wave input; that is, the frequency 
components of the sounds are characteristic of the measuring system 
and are not characteristic of the dynamics of the blood-filled 
artery. This does not mean that the sounds produced are independent 
of the properties of the blood-filled artery, but the dependence is 
only on the sharpness of the pulse (which depends on the properties 
listed in the mechanism hypothesis) which excites the measuring 
system. 

The term "measuring system” means the ear alone or the aided 
ear (i.e., with a stethoscope, or some microphone speaker system). 

Sounds, by definition, require the ear as part of the sound 
measuring system. 

This hypothesis implies that the "frequencies" of the blood artery system 
cannot be obtained directly from the frequency spectral analysis of the sounds. 
In fact, those frequencies obtained from such- analyses reflect the inherent 



frequencies of the measuring system. This could account for the vide 
variation of "frequencies " for the Korotkoff sounds reported in the literature. 


Experimental Corroboration 

Although Erlanger and Bramwell had meager analytical support of their 
hypothesis, they did extensive experimental work. A summary of their work on 
mechanical (as opposed to biological) systems was given by Bramwell (ref. 32) 
as follows: 

" Sounds produced by waves in closed tubes . How is it that the 
deflation of the armlet causes the changes in the arterial sounds 
described by Korotkoff? To answer that question I must first refer 
to some experimental observations made on an artificial schema of 
the circulation. This consisted of the inner tube of a bicycle tyre 
containing fluid but no air. The distal end of the tube was closed 
by a wooden bung and the proximal end connected with a pump capable 
of generating waves in the fluid. The tyre represented a portion of 
the arterial tree and the pump the heart. Erlanger (1924) found 
that a wave passing along a rubber tube filled with water produced 
a sound which varied both in quality and in intensity at different 
points along the tube. The change in quality of the sounds, as 
heard with a stethoscope, enables one to map out three regions in 
the tube. 

(1) In the proximal region, only a dull thud is heard; this 
becomes louder as one passes farther down the tube. 

(2) In the intermediate region, the sound is very loud; here it 
is no longer dull but sharp like the crack of a whip. 

( 3 ) In the distal region, the sound loses its cracking character 
and becomes progressively fainter as we pass towards the distal end 
of the tube. 

"Erlanger pointed out that, in the intermediate region of the 
tube where the loud crack is heard, the wall of the tube is subject 
to a severe shock. This he demonstrated by sprinkling fine sand on 
the upper surface of the tube. With the passage of a wave the sand 
is violently thrown off from that region in the tube where the shock 
is most severe and where the loud cracking is heard on auscultation, 
but lies undisturbed in the segments on either side. 

"Working in Erlanger f s laboratory in 1923; I made some further 
observations on this phenomenon (Bramwell 1937)* I found that the 
position of the disturbed area - i.e., the area in which the sand is 
thrown off the tube - varied with the pressure of the fluid in the 
tube. As one raises the pressure, the disturbed area recedes from 
the proximal towards the distal end of the tube, at high pressures 
the wave passes down the whole length of the tube without producing 
sufficient disturbance to throw off the sand. In other words, the 
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higher the initial pres sure , the farther has the wave to travel 
before the conditions arise which lead to the throwing off of sand 
and the production of sharp cracking sounds." 

Bramwell presented sufficient detail of his experiments (ref. 20) for a 
comparison with the theory of section II. Using mercury as the fluid, 

Bramwell introduced the pressure wave shown in figure 15. He determined 
directly the wave velocity pressure curve from the distortion of the wave. 

The result is shown in figure 1 6 . 

The calculation of the critical length of tube (eq.. (56)) is particularly 
simple for the given pulse wave (almost linear in its rise) and the minimum 
critical length is also easily obtained (at minimum c 2 /(dc/dp)). If it is 
now assumed that the minimum critical length corresponds to the "disturbed 
area" determined experimentally by Bramwell, a comparison of theory and experi- 
ment is in order. Because of apparatus restrictions, Bramwell could determine 
only bounds on the "disturbed area" and these are shown as vertical lines in 
figure 17* The calculated minimum critical length is shown as a solid line. 



Figure 15 •- Experimental input pulse wave (Bramwell data). The pressure is the differential 

pressure across the tube 'wall. 


^3 




Figure 17*- Comparison of experimental (Branrwell) and theoretical (this report) determination of 
critical tube length. The pressure is the differential pressure across the tube vail. 


V. RESUME 


The results of the analytical and experimental work presented in this 
report can “best he summarized by a series of conclusive statements. 

1. A solution for the finite amplitude -wave equation with general tube 
pressure-area relation can be obtained for simple (nonreflecting) waves. 

2. The finite amplitude nonlinear wave equation reduces to the linear 
wave equation when the pressure -area relation for a longitudinally constrained 
rubber tube with circular cross section is introduced. 

3. The set of criteria for wave distortion (steepening or nonsteepening) 
can be stated as a function of the pressure-area relation of the tube. 
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4. The distortion of vaves and the corresponding critical length of tube 
can be calculated if the initial shape of the wave (entering the tube) and the 
pressure-area (or wave velocity-pressure) relation of the tube are known. 

5- The development of shock waves In fluid -filled elastic tubes is theo- 
retically feasible (from the simplified theory). The shock -wave velocity and 
energy loss are readily obtainable if the amplitude of the wave and the 
pressure -area relation for the tube are known. 

6 . The analysis in this report supports the wave distortion hypothesis 
and experimental work of Erlanger and Bramwell. 

Statements 1 } 8 , and 9 are based on new hypotheses presented in this report. 
These hypotheses have not been proven but are supported by the analysis and 
experiments discussed in this report. 

7. The development of steep waves produces sharp waves or possibly shock 
waves and not "breakers” as hypothesized by Bramwell. 

8 . In the indirect method of blood pressure measurement ( sphygmomanome- 
try) the cuff pressures at which Korotkoff sounds are heard correspond to 
artery conditions suitable for the pulse wave to steepen sufficiently to 
excite a response in the measuring system. 

9 . The frequency components of the Korotkoff sounds are characteristic 
of the measuring system (stethoscope or other) and not characteristic of the 
blood-filled artery. More precisely, the pressure variation inside the artery 
is very nearly a step (with possibly a "blip" or two) while the pressure varia- 
tions to which the ear is exposed (through a stethoscope or otherwise) contain 
oscillatory components of the measuring system response (to the steep wave). 


Ames Research Center 

Rational Aeronautics and Space Administration 
Moffett Field, Calif., 94035, May 27, 1968 
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